Doping of graphene by a Au(lll) substrate: Calculation strategy within the local 
density approximation and a semiempirical van der Waals approach 
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We have performed a density functional study of graphene adsorbed on Au(lll) surface using both 
a local density approximation and a semiempirical van der Waals approach proposed by Grimme, 
known as the DFT-D2 method. Graphene physisorbed on metal has the linear dispersion preserved 
in the band-structure, but the Fermi level of the system is shifted with respect to the conical points 
which results in a doping effect. We show that the type and amount of doping depends not only on 
the choice of the exchange-correlation functional used in the calculations, but also on the supercell 
geometry that models the physical system. We analyzed how the factors such as the in-plane cell 
parameter and interlayer spacing in gold influence the Fermi level shift and we found that even 
a small variation in these parameters may cause a transition from p-type to n-type doping. We 
have selected a reasonable set of model parameters and obtained that graphene is either undoped 
or at most slightly p-type doped on the clean Au(lll) surface, which seems to be in line with 
experimental findings. On the other hand, modifications of the substrate lattice may induce larger 
doping up to 0.30-0.40 eV depending on the graphene-metal adsorption distance. The sensitivity 
of the graphene-gold interface to the structural parameters may allow to tune doping across the 
samples which could lead to possible applications in graphene-based electronic devices. We believe 
that the present remarks can be also useful for other studies based on the periodic DFT. 

PACS numbers: 73.22.Pr, 71. 15. Mb, 73.20.Hb 
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I. INTRODUCTION 

Graphene, a two-dimensional sheet of carbon atoms 
packed into a honeycomb lattice, has attracted con- 
siderable attention due to its thin geometry and un- 
usual electronic properties promising for a wide range 
of applications One of the most important themes in 
recent research is to investigate the interaction between 
graphene and its surrounding environment, which is of 
practical relevance for the process of samples production 
and fabrication of graphene-based electronic devices. In 
particular, metallic substrates are used as catalysts in 
graphene formation, as probes during the electrical mea- 
surements and as source/drain electrodes in electronic 
devices. 

Early systematic density functional theoretical stud- 
ies based on the local density approximation (LDA) have 
shown that two types of interfaces can be formed between 
graphene and metal, i.e. physisorption and chemisorp- 
tion should be distinguished* 4 -! 5 - The bonding of graphene 
to Al, Ag, Cu, Au and Pt(lll) surfaces is weak and pre- 
serves graphene's characteristic Dirac cones, while Co, 
Ni and Pd bind graphene strongly which disturbs its ir 
bands through coupling with the d-orbitals of the met- 
als. The physisorption of graphene on the metal surface 
causes a Fermi level shift downward (upward): It means 
that holes (electrons) are donated by the metal substrate 
to graphene which becomes p-type (n-type) doped. Ac- 
cording to this theoretical study, graphene on Al, Ag and 
Cu is n-type doped, while the interaction with Au and Pt 
(111) surfaces causes p-type doping. Most of these LDA 
results were confirmed in experiments^— 

On the other hand, the local density approximation is 



well known to overestimate the binding in weakly bonded 
systems where interactions are mainly of van der Waals 
(vdW) type. It is regarded to be an unreliable method 
for inhomogeneous systems, such as graphite and many 
organic compounds 1,2 . A generalized gradient approxi- 
mation (GGA) does not improve the results. It tends 
to underestimate the binding: in the case of graphene 
on gold it gives no binding at all, while for graphene on 
Ni only a metastable minimum is predicted. Although 
GGA in general offers reasonable predictions for met- 
als, it gives incorrect binding behavior for the graphene- 
meta l 12 ' 13 and for other organic-inorganic interface s 14 ' 15 . 
Thus, much research effort has been devoted to include 
vdW forces in the studies of graphene-metal composites. 

In the paper of Vanin et al.— the recently developed 
van der Waals density functional (vdW-DF) 16 ' 17 has been 
used to study the metal-graphene interfaces. Contrary 
to the previous LDA results, the inclusion of non-local 
correlations leads to weak binding for all metals and 
gives graphene-substrate separations in the range of 3.40- 
3.72 A. At these distances graphene is physisorbed and 
the graphene band structure is unaffected by the sub- 
strate: Strong binding in the case of Ni, Co and Pd is 
not predicted by the vdW-DF, thus there is no band-gap 
opening at the K-point. These results are in a strong 
conflict with low-energy electron diffraction (LEED) 18 
and angle-resolved photoemission spectroscopy (ARPES) 
measurements 7 - for graphene on Ni, which predict an ad- 
sorption distance of 2.11 A and the presence of a band- 
gap. It seems that the DFT calculations should be tested 
first for the graphene-Ni interface so that the adsorption 
distance and band structure are correct. 

Furthermore, in a very recent study 19 of graphene ad- 



2 



sorbed on (111) surfaces of Ni, Cu, Pd, Ag, Au and Pt 
the second version of vdW-DF (vdW-DF2)2° with a C09 
exchange functional developed by Cooper 21 has been ap- 
plied and assessed to be the most appropriate combina- 
tion of exchange-correlation interaction. It gives accu- 
rate adsorption geometries and electronic structures in 
accordance with available experimental data for most of 
the metals. Two classes of interfaces, physisorption and 
chemisorption, are distinguished, in agreement with pre- 
vious LDA results. The band-gap opening at the K point 
observed on the Ni (111) is reproduced reasonably well. 
The authors also compare their graphene-copper distance 
with semiempirical DFT-D (Ref|22j) calculations and con- 
clude that both methods give nearly the same equilibrium 
distances. 

One important discrepancy when comparing to previ- 
ous theoretical and experimental studies, is reported in 
Reffl9l. The obtained n-type doping of graphene on gold 
contradicts the experimental findings of Klusek et al., 
who observed that holes are donated to graphene on Au 
substrate. The p-type doping has been confirmed in the 
high-resolution ARPES measurements of Varykhalov et 
al££. It has been demonstrated that graphene on Ni in- 
tercalated by a Au monolayer remains gapless and that 
a very small hole doping (AEp — +100 ± 20 meV) is 
observed in the band structure. Moreover, it has been 
shown that also the presence of gold atoms on graphene 
causes its p-type doping^!. In the Reffl9l. it is suggested 
that the origin of this inconsistency for the graphene- 
Au system is connected with the sensitivity of the dop- 
ing polarity to the adsorption distance. It seems that a 
more accurate description of the graphene-gold contact 
is needed. 

Theoretical research concerning graphene on the met- 
als' surfaces has been focused on the construction and 
choice of the appropriate exchange-correlation function- 
al. They should provide a correct determination of 
the graphene-substrate separation and, as a consequence, 
capture the main characteristic features of the spectrum, 
as an energy gap in the band structure of chemisorbed 
graphene as well as the type and amount of doping for 
weak adsorption. The correct determination of the ad- 
sorption distance is especially crucial, when the transi- 
tion between extremely different properties occurs near 
the estimated distance, in particular the change from n- 
type to p-type doping of graphene on gold 5 . 

However, there are other factors that are important 
in the modeling of interfaces, especially those related to 
the geometry of the supercell. First, performing peri- 
odic DFT calculations requires direct matching of the 
graphene's and substrate's unit cells. Typically, the C- 
C bonds of physisorbed graphene are not stretched or 
quenched in the physical sense. On the other hand, ad- 
justing the substrate to the graphene generates strain in 
its lattice, because the interatomic distances in a metal 
are artificially contracted. It results in a change in the 
surface relaxation and can additionally alter the elec- 
tronic properties. In the case of a slight lattice mismatch, 



these two competing geometries are considered equiva- 
lently in literature 5 -^. Furthermore, there are various 
routes for the structural optimization of the graphene- 
metal system. Relaxation leads to different geometries 
and electronic properties if i) the metal atoms are kept 
in their experimental posistions 12 , ii) only carbon atoms 
are fixed 2 - 5 -, or iii) graphene and the two top layers of gold 
are allowed to move 5 -. The above-mentioned factors give 
small percentage changes in results for almost all types of 
the substrates. However, the graphene-gold interface is 
very sensitive to variations in the adsorption parameters 
and qualitatively different properties can be derived. 

In this paper, we attempt to resolve the issue of doping 
graphene on Au(lll). In our DFT calculations both lo- 
cal density and generalized gradient approximations have 
been used. To include vdW corrections, we have ap- 
plied the semi-empirical DFT-D2 Grimme's metho d 22 ' 26 , 
which has brought satisfactory predictions for similar sys- 
tems. For example, the adsorption of aromatic molecules 
on the (111) surfaces of noble metals^ and the adsorption 
of Cu, Ag and Au atoms on graphene 2 ^ were accurately 
described. DFT-D2 should be then a valuable tool also 
for graphene-metal contact studies. 

First, we have performed a series of fixed geometry 
calculations to clearly identify the main model factors 
determining the doping of graphene on gold. We have 
calculated the total energies as well as the shifts of the 
Fermi level with respect to the Dirac points using two dif- 
ferent functionals. We show that the choice of exchange- 
correlation interaction approximation alters the doping 
directly and also indirectly by influencing the adsorption 
distance. Second, we have evaluated the two-dimensional 
maps of Fermi level shifts function of graphene- 

gold separations and the in-plane unit cell parameters 
assumed for the simulations in order to identify regions 
where graphene is n-type doped, p-type doped or un- 
doped. We have analyzed how the Fermi level shift varies 
with a change of interlayer spacings in gold, which have 
also significant physical implications. Finally, we have 
calculated surface relaxation and studied its influence on 
electronic properties. This detailed analysis allows to as- 
sign the inconsistency indicated in Ref. [TH to geometry 
parameters defined in the model. 

The main purpose of this paper is to develop a method- 
ology, so that the calculations with vdW corrections give 
correct results for all typical metals and to propose a 
set of model parameters which are reasonable from both 
computational and experimental points of view. Since 
the recommended choice of interface geometry prevents 
a reliable relaxation of the Au(lll) surface, alternative 
routes to manage this limitation are suggested. More- 
over, in order to demonstrate the reliability of a simple 
DFT-D2 method, we have performed basic calculations 
of a graphene-nickel interface. 

The present paper is organized as follows. In Sec. II 
the calculation method and computational details are 
summarized. In Sec. Ill we present the interpretation 
of fixed-geometry calculations and identify factors deter- 
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FIG. 1: (Color online) Side view (a) and top view (b) of 
adsorption geometry of graphene on Au(lll). Carbon atoms 
are denoted as black balls, gold atoms as yellow (gray) ones. 
Parallelogram defines the unit cell. 

mining the doping. Surface relaxation is analyzed in Sec. 
IV, while conclusions and perspectives are discussed in 
Sec. V. 



II. METHODS 

We have performed ab initio DFT calculations us- 
ing the VASP cod o 31 ' 32 equipped with the projector 
augmented wave (PAW) metho d 33 ' 34 for electron-ion 
interactions. The exchange-correlation interaction is 
treated in the generalized gradient approximation in 
the parametrization of Perdew, Burke and Ernzer- 
hof (PBE)2£. Long-range dispersion corrections have 
been taken into account within a DFT-D2 approach 
of Grimme 2 ^, as implemented in the latest version of 
VASP 26 . Since a dispersion coefficient C§ for gold has 
not been listed in the original paper of Grimme 2 ^, we 
have used a value of 40.62 Jnm 6 /mol and of 1.772 A for 
the vdW radius of Au (-Ro), according to the suggestion 
given in Ref. [2^. Also, as in Ref. [28|, the pair interactions 
up to a radius of 12 A have been included in the calcula- 
tions and the global scaling factor sq has been set to 0.75 
because the PBE functional was chosen. For compari- 
son, most of the calculations have been also done using 
the local density approximation (LDA)2&. The electronic 
wave functions have been expanded in a plane-wave ba- 
sis set of 400 eV. The electronic self-consistency criterion 
has been set to 10~ 7 eV. 

The Au(lll) surface is modeled by a periodic slab ge- 
ometry: each supercell contains six atomic layers of metal 
and a graphene sheet adsorbed on one side. A vacuum 
spacing of at least 13 A is used in the direction normal to 
the interface in order to avoid interactions with spurious 
replicas. We employed the in-plane adsorption geometry 
suggested in Refs.00 which is illustrated in Figjl] where 
2x2 graphene and substrate (yo x Vo) unit cells are di- 
rectly matched. We have performed calculations for both 
configurations, i.e. with the lattice constant of the metal 



adapted to that of graphene and with the graphene lat- 
tice constant adjusted to the substrate one. The in-plane 
cell parameter a n ex can be defined as: 

G o G Au Au % 

"hex = 2 a . "hex = « Y 2 W 

where a G and et Au stand for the lattice constant of 
graphene and bulk gold, respectively. 

Consistently, the Au(lll) interplane distances can be 
calculated from: 



However, this can give too short distances if the metal is 
adjusted to the graphene, and it seems to be more reason- 
able to set it independently, especially when relaxation is 
not included. The above-mentioned parameters as well 
as their relation to the lattice constants are illustrated 
in Fig[TJ Note that in the case of structural optimiza- 
tion the parameter <jj nt is no longer a constant. Instead, 
there appear different values d\ 2 , g?23, etc. We define d 
as the fixed graphene-substrate distance and dc-Au as 
optimized distance between graphene and gold. 

Moreover, graphene's and metal's lattice constants can 
be determined according to various approaches. They 
can be either optimized using LDA and PBE functionals, 
or measured experimentally. It gives a wide choice of 
values that determine the size of the supercell (ahex) and 
the corresponding interlayer distances d lrA . Typical sets 
of these parameters calculated from different values of 
a G and a Au are listed in Table I. We have performed 
calculations for most of these sets. 

The DFT-D2 method allows to include vdW correc- 
tions within each ionic step. During the structural opti- 
mizations done with a conjugate gradient algorithm the 
top two layers of gold as well as all carbon atoms have 
been relaxed. Total energies were converged to within 
10~ 6 with respect to the ionic steps. 

For accurate Brillouin zone integrations we use the 
tetrahedron scheme 3 -^ and the T-centered 24x24x1 k- 
point mesh. 

As a benchmark test of the PBE+D2 approach applied 
to graphene/metal interfaces, we first optimized the ge- 
ometry of graphene/Ni system and determined its elec- 
tronic properties. All parameters needed in calculations 
are the same as those described above. (The Cq and Ro 
coefficients for Ni arc defined in VASP.) Moreover, the 
spin polarization is taken into account. The supercell 
geometry is set as proposed in ReflU, where the nickel 
(1 x 1) unit cell is adjusted to graphene with a LDA- 
optimized lattice constant equal to a = 2.445 A. The 
LDA calculations of the structure lead to the adsorption 
distance of 2.05 A , thus the results given in RefJa are 
well reproduced. In contrast, the PBE+D2-optimized 
graphene-nickel distance achieves nearly 2.15 A (2.11 
without spin polarization) that is still in line with exper- 
imental findings 1 ^. A band-gap at the graphene K-point 
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TABLE I: The sets of parameters determining the size of the supercell and the geometry of the slab. For bulk Au the 
LDA-optimized lattice constant compare well with the experimental valued. The one obtained using the PBE functional is 
overestimated, but in accordance with previous DFT results^. The interlayer spacings in gold dj nt are calculated from Eq. (2). 
All values are given in A. 
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is comparable to the one obtained within local density 
approximation, which is smaller than the one observed 
in the ARPES experiment 7 . The use of vdW-DF2 C09 
might give more accurate values, but we believe that the 
DFT-D2 method is sufficiently precise for a description 
of graphene-metal systems. 



III. FIXED GEOMETRY CALCULATIONS: 
MAIN FACTORS DETERMINING THE DOPING 

The weak binding of graphene on gold preserves the 
characteristic conical ir bands touching at the K/K' 
points. In contrast to the free-standing layer, the Fermi 
level is shifted with respect to the Dirac points due to the 
interaction with the metal surface. We define the doping 
level as AEp — Ed — Ep, where Ed is the energy where 
the Dirac cones touches and Ep denotes a Fermi energy 
of the system. In FigJ5] we present the band-structure 
of free-standing graphene and of graphene on gold calcu- 
lated in a fixed geometry using a local density approxi- 
mation. In the latter case, the Fermi level shift AEp is 
well visible and estimated to be about 0.20 eV. The band- 
structure is very similar to the one presented in Fig. 2 of 
RefH in spite of the fact that relaxation has not been 
included. Overall, the whole dispersion hardly changes 
with the aforementioned model parameters: only a dop- 
ing level AEp can be altered significantly. We have then 
performed fixed-geometry calculations, since it allows to 
discuss the role of each parameter separately. 



A. Graphene-substrate distance 

It has been previously noted that the doping level and 
polarity are very sensitive to adsorption distance. A vari- 
ation in the Fermi level shift as a function of graphene- 
substrate separation calculated within a LDA is shown 
in Fig[3] (black upper line). Similar dependencies have 
been presented in Fig. 6 of Ref|5| for all typical metallic 
substrates. One can observe^ that graphene on Al(lll) 
is estimated to be n-type doped for a whole range of 
distances, while for the Pt substrate simulations consis- 
tently yield p-type doping. On the other hand, for Ag, 
Au and Cu the change in doping polarity from n-type 
to p-type occurs in the vicinity of the estimated adsorp- 
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FIG. 2: (Color online) Top panel: The electronic structure of 
free-standing graphene obtained with a 2 x 2 unit cell. The 
inset shows a Brillouin zone with the high-symmetry points 
labeled by F, K and M. Bottom panel: The band structure 
of graphene on Au(lll). The red lines denote the graphene 
bands with a n character. The zero energy is at a Fermi level, 
d int =2.35 A, a = 2.445 A, d = 3.31 A. 
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FIG. 3: Fermi level shift with respect to the Dirac point 
as a function of graphene-gold distance d. The values of 
Qhcx = 4.89 A and d lnt = 2.35 A have been set for the cal- 
culations. The black (upper) curve has been obtained with a 
LDA approximation and the red one with a PBE functional. 
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FIG. 4: Binding energy (Eb) per carbon atom of graphene on 
Au(lll) as a function of graphene-gold distance d calculated 
with a LDA (upper, black line) and the PBE+D2 method (red 
line). The parameters are set to ahex=4.89 A and d = 2.35 

A. 



tion distances. Moreover, as one can easily see in Fig|3] 
and in Ref||| (Fig. 6), such dependencies vary signifi- 
cantly in this region, i.e. in the vicinity of the transition 
the slope of the curve is significant, which explains an 
extreme sensitivity of doping to the graphene-metal sep- 
aration. Both effects may inhibit reliable determination 
of the type and level of doping, especially if the distance 
has not been measured experimentally and its estimation 
relies only on DFT calculations. 

We have evaluated the same dependency using the 
PBE parametrization for a description of the exchange- 



correlation interaction (the red bottom line in Fig[3]). 
The vdW corrections are added to the conventional 
Kohn-Sham DFT energies, thus they do not alter the 
band-structure^. It should be stressed that the doping 
level differs from that obtained using the LDA, even if 
the issue of adsorption distance is ignored. The lines 
in Fig[3] illustrate that there are regions (3.15-3.23 A), 
where the LDA and PBE provide opposite doping po- 
larities, although the numerical difference is small. A 
similar comparison has been presented in Fig. 9 of Ref.5 
for a graphene-copper system. However, the Fermi level 
shifts AEp calculated with the GGA were within 0.07 
eV of the ones obtained with the LDA, which is irrelevant 
if it does not lead to qualitative changes. 

Obviously, the adsorption distance depends on the 
choice of the functional. It can be easily noticed in 
FiglU where the binding energies as a function of distance 
have been plotted. The black (top-most) line denotes the 
LDA-obtained results, whereas the PBE+D2 predictions 
are given by the red curve. The PBE+D2 calculations 
lead to an adsorption distance of 3.22 A, while the LDA 
estimation is d — 3.36 A. This gives a difference in doping 
from zero to more than +0.25 eV (see Fig[3]), and thus 
both direct and indirect influences of the functional lead 
to a significant general change. 



B. In-plane lattice constant 

In Sec. II we have indicated that the size of the super- 
cell parametrized by ahex can be set to any value from 
those listed in Table I. In literature, the extreme values, 
4.89 and 5.11 A, have been used equivalently. However, 
the physical meaning of this parameter is a stretch of 
graphene bonds, although they should be rather unaf- 
fected by weak adsorption. 

We have calculated dependencies similar to those pre- 
sented in Fig J3] for a few values of ahex and interpolated 
them into a whole range of in-plane lattice constants 
(4.89-5.11 A). In FiglUwe present the LDA-obtained two- 
dimensional map of Fermi level shifts AEp(d, ahex)- The 
color represents a value of doping for a given pair of in- 
plane cell parameter and graphene-gold distance (ahex, 
d). One can easily see that graphene could be n-type 
doped if the C-C bonds were stretched sufficiently It is 
clear, however, that in the case of LDA calculations set- 
ting such values of ah cx is, in fact, unreasonable, because 
LDA optimized values are ahex=4.89 A or ahcx=4.976 A 
(see Table I). It is worthwhile to notice that even tak- 
ing the experimental ah cx =4.997 A value, results in a 
considerable decrease in the value of p-type doping. 

The influence of these effects has been carefully stud- 
ied in Refi for copper substrate due to the fact that 
the lattice mismatch between graphene and the Cu(lll) 
surface is the largest among the substrates considered in 
this study. A comparison between configurations with 
copper's lattice constant adapted to graphene and with 
graphene stretched to the in-plane lattice constant of Cu 
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FIG. 5: (Color online) The Fermi level shift relative to the 
conical point as a function of in-plane lattice constant «hex 
and graphene-gold distance d. The doping changes from n- 
type (darker) to p-type (lighter). The white line denotes a 
contour of zero doping. The right-hand panel shows the color 
scheme used in the map. Calculations has been done using the 
LDA approximation and spacing di nt =2.35 A (a) and using 
the PBE parametrization and interplane distance di n t=2.35 

A (b) 



gives errors within 0.15 eV. It seems that this result could 
influence only the value of n-type doping, because the 
transition of polarities occurs for distances above the op- 
timized one in the case of both configurations (see Fig. 9 
in Ref.5). However, for the gold substrate a relative shift 
of 0.15 eV can lead to more considerable changes or even 
cause a doping type transition. It seems that the lat- 
tice mismatch should not be the only criterion for such 
investigations. 

A similar map calculated with the PBE functional is 
shown in FigfS] (b). It can be noticed that the variation 
in doping is even stronger in a relevant range of d, which 
is related to the fact that the PBE dependency lies be- 
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FIG. 6: (Color online) The Fermi level shift relative to the 
conical point as a function of graphene-metal distance for 
Qhcx = 4.89 A and for different values of d lnt : the blue (bot- 
tom) line is assigned to di nt = 2.305, a red one (next) to 2.345 
A, the green one to 2.41 A and the purple (the top-most) one 
to 2.87 A. The PBE functional has been used in the calcula- 
tions. 

low the LDA one (Fig[3]). We would like to stress that for 
ahcx=5.11 A, as chosen in Ref.19, the electrons are always 
donated to graphene from the Au substrate resulting in 
n-type doping. We have also performed similar calcula- 
tions with a slightly increased value of d; nt . It yields a 
systematic shift of doping, but does not invalidate the 
overall conclusion. 



C. Interlayer spacing in gold 

According to the remarks given in Sec. II, we can 
choose the interlayer spacing in gold <ii nt from a range of 
values defined in Table I. The choice of the starting dj nt 
is less important if the structure is allowed to relax. Since 
locking the in-plane lattice constant may lead to incor- 
rect optimization, we assume that it is better to study 
first the configurations with fixed di n t spacings chosen 
around the experimental bulk values. 

Figure [S] illustrates the dependencies AEp (d) for dif- 
ferent values of d- m t- The lines in the middle represent 
typical choices for the interplane spacing (see Table I). 
We conclude that the larger values of d- m t lead to an 
increase in the shift of the Fermi level. 

The most interesting feature of the plot in Figgis the 
top-most (purple) line representing the function AEp(d) 
for a large value of c?i nt . First, it allows to clearly iden- 
tify the direction of the doping changes with an increase 
in the value of dj n t. The properties of such an artificial 
surface also suggest how strongly the substrate modifi- 
cations may affect the doping. It seems in this case that 
for distances d higher than 3.4 A the doping saturates 
and achieves a value of AEp=0A0 eV. It means that any 
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crystal imperfections of the substrate can induce changes 
in doping from zero to about 0.30-0.40 eV depending on 
the adsorption distance. 



IV. SURFACE RELAXATION 

DFT studies of graphene/metal interfaces typically in- 
clude structural optimization. The surface relaxation not 
only could alter the doping level, but also is crucial for the 
evaluation of phonon dispersion and (scanning tunneling 
microscopy) STM topographies. In the previous section, 
we have reported the results of fixed-geometry calcula- 
tions in order to analyze each model parameter separately 
and study its influence on the electronic properties. Now, 
we show that the relaxation of compressed unit cell of the 
substrate could lead to an unphysical expansion of its in- 
terlayer spacings. It means that performing calculations 
in the configuration with at least the metal atoms fixed 
might give more reasonable results. 

The structural optimization of a system consisting of 
graphene and metal can be done by many ways. Ac- 
cording to suggestion given in Refi we relax all carbon 
atoms and the top two layers of gold until the total energy 
change between subsequent steps is smaller than 10~ 6 
cV. It should be mentioned that other approaches are 
frequently used: the relaxation of carbon atoms only 12 , 
relaxation of gold atoms only 2 ^, relaxation of graphene 
and gold separately followed by distance optimization 19 
etc. Obviously, none of these methods can take into ac- 
count the unique herringbone reconstruction of the clean 
(111) surface of Au22r— . We have demonstrated that 
in-plane optimization hardly affects the band-structure, 
however the influence of reconstruction is very difficult to 
estimate within DFT with periodic boundary conditions 
and was not taken into account. 

We focus on the changes in the interlayer spacing in 
gold upon relaxation, since we have shown that it could 
influence the doping level AEp. We study the relax- 
ation effects in metal taking into account first d\2 and 
second (^3 interplane distance in gold. We have per- 
formed structural calculations using both LDA and PBE 
exchange-correlation functionals. The results for various 
configurations of starting geometry are listed in Table II. 
In particular, the relaxation of a clean Au(IlI) surface, 
calculations with a larger unit cell, as well as optimization 
of only carbon atoms have been performed (see the cap- 
tion of Table II). First, it should be stressed that match- 
ing gold to graphene leads to considerable expansions of 
d±2 and c?23, even for clean Au(lll) surface. Although 
the LDA and PBE+D2 methods give quantitative differ- 
ences, for both functionals the values of interlayer spac- 
ings d\2 and (^23 decrease after setting greater values of 
ahcx- The interlayer distances calculated for ahox=5.02 A 
are close to the ones obtained in the previous work based 
on the GGA approximation— and seem to be the most 
reasonable values. On the other hand, this configura- 
tion proves to be inconsistent with experimentali - n-type 



doping due to the stretching of C-C bonds, as discussed 
in the previous section. It is clear that the choice of 
«hcx which allows to preserve the length of the graphene 
bonds prevents reliable structural optimization. The ver- 
tical expansion of gold causes also an increase in doping, 
which could be equally extrapolated from data in FiglHl 
In our opinion, the ahex parameter adapted to graphene 
needs to be preserved to ensure correct electronic proper- 
ties of graphene, but the expansion of interlayer spacings 
in metal should be prevented by keeping its atoms in their 
positions previously optimized with larger cthex- In fact, 
limiting the degrees of freedom of the DFT calculation 
by fixing the in-plane lattice constant can lead to errors 
comparable even to those caused by locking the whole 
system. 



V. DISCUSSION AND PERSPECTIVES 

The theoretical study reported in the previous sections 
has been focused on the analysis of the model depen- 
dent factors that strongly affect the obtained physical 
properties of graphene on gold. Indeed, the in-plane lat- 
tice constant has been treated as a free parameter, al- 
though its physical meaning is related to the length of the 
graphene's bonds. However, we estimated that the ad- 
sorption distance hardly depends on ahcx suggesting that 
the interaction with the gold substrate should not con- 
siderably stretch the bonds. On the other hand, available 
experimental data indicate p-type doping of graphene on 
a gold substrate 10 . According to our studies, it means 
that the C-C bonds do not change their length. In or- 
der to obtain correct electronic properties of graphene on 
Au(III), the lattice constant of graphene should rather 
be preserved. Effectively, it represents the case where 
metal atoms with a given coverage were deposited on 
graphene^. 

It should be noted that in most of our PBE+D2 cal- 
culations we have set ahex=4.89 A consistently with the 
LDA-optimized lattice constant of graphene, just to eas- 
ily compare the results with previous works based on 
LDA. According to the data presented in Figj5](b), set- 
ting a real experimental graphene lattice constant (a = 
2.46 A) should provide very slight n-type doping. How- 
ever, the simulations with ahcx=4.932 A lead to nearly 
zero doping due to the compensation by relaxation ef- 
fects. Overall, calculations performed within the DFT- 
D2 method indicate that graphene is undoped on a clean 
perfect Au(lll) surface meaning that it would be nearly 
neutral (note that DFT-D2 tends to underestimate the 
adsorption distance) . It could be an explanation why the 
intercalation of graphene-Ni system by gold atoms leads 
to almost free-standing graphene. 

We have analyzed the graphene-gold system because 
contradictory data have been reported in the litera- 
ture and we demonstrated that these inconsistencies are 
hardly connected with the choice of the functional. The 
case of a gold substrate can appear to be complicated, 



TABLE II: Optimized structural parameters and the calculated Fermi level shifts AEf in graphene-gold system for different 
starting configurations (cthex, di nt ). In the first column with the slashes, we present also the results for a clean Au surface. 
The fourth column with the slashes, beside standard data, contains also the results of relaxation of only carbon atoms and of 
a clean metal surface, respectively. All distances are given in A. 



LDA PBE+D2 

a hex 4.89 4.976 4.89 4.89 5.02 4.89 

dint 2.345 2.345 2.41 2.345 2.345 2.41 

dc-Au 3.365/- - 3.35 3.21/3.22/- 3.20 3.20 

di2 2.43/2.425 2.36 2.425 2.52/2.345/2.55 2.43 2.54 

d 2 z 2.39/2.39 2.34 2.39 2.44/2.345/2.45 2.36 2.46 

AE F (eV) 0.26/- - 0.25 0.08/0.00/- -0.16 0.07 



because the transition of n-type to p-type doping is 
abrupt and occurs in the region of the estimated ad- 
sorption distance. However, it seems that the structure 
and electronic properties of graphene on other metallic 
substrates, as for example on copper, can be more com- 
plex, because the adsorption distance has been reported 
to change with ahexr^- We have used the example of 
gold to unravel the importance of structural details re- 
sponsible for changes in the Fermi level shift. It is clear 
that similar methods might be used for characterization 
of graphene on other metals, especially on Cu and Ag. 

Moreover, our conclusions could be also useful for 
further studies of molecular layers on the noble metal 
substrates. Gold surfaces are frequently employed 
as templates for physisorbed self-assembled monolayers 
(SAMs)2ir— . Since various choices of geometry and 
structural optimization methodologies are applie d 45 ' 46 , 
our results indicate that similar problems can occur in 
these systems, especially for the determination of the 
molecule (assembly)-substrate distance, the values of 
charge transfer and binding energies. 

Finally, we have checked the validity of DFT-D2 
method applied to graphene/metal interfaces. We have 
studied graphene chemisorbed on Ni as well as ph- 
ysisorbed on a Au surface. The results are more 
reasonable than those obtained within the vdW-DF 
functional 12 , especially the estimation of a graphene- 
substrate distance is more satisfactory. However, the 



binding energies seem to be considerably overestimated, 
which was also reported in studies of the molecules on 
metals surface s 15 ' 27 . Of course, a systematic study is re- 
quired to fully confirm the accuracy of this method. 

We believe that the use of the vdW-DF C09 functional 
could be equally efficient, since the results in conflict 
with the experiment reported in Refill might have been 
caused by the choice of the in-plane lattice constant (see 
FigJS]). This configuration seems to describe the samples 
different from those studied in experiments. 
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